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ABSTRACT 

We present a new determination of the co-moving star formation den¬ 
sity at redshifts z ^ 0.35 from the 1.4 GHz luminosity function of sub-mJy 
star-burst galaxies. Our sample, taken from Benn et al. (1993), is insen¬ 
sitive to dust obscuration. The shape of the Luminosity function of this 
sample is indistinguishable from a number of reasonable a prior models 
of the luminosity function. Using these shapes we calculate the modest 
corrections (typically ^ 20 per cent) to the observed 1.4 GHz luminosity 
density. We find that the cosmic variance in our estimate of this lumi¬ 
nosity density is large. We find a luminosity density in broad agreement 
with that from the RSA sample by Condon et al. (1987). We infer a co¬ 
moving star formation rate surprisingly similar to coeval estimates from 
the Canada-France Redshift Survey, in both ultraviolet and Ha, although 
the later may also be affected by Cosmic variance. We conclude that the 
intermediate 0.05 < z < 0.3 star-formation rate is not yet well determined 
partially due to uncertain extinction corrections, but also partially due to 
cosmic variance. We suggest that deep moderate area radio surveys will 
improve this situation considerably. 

Key words: galaxies: formation - infrared: galaxies - surveys - galaxies: 
evolution - galaxies: star-burst - galaxies: Seyfert 


1 INTRODUCTION 

In recent years a major focus of extra-galactic astro¬ 
physics has been in estimating the evolving, volume- 
averaged star formation rate. This work has been par¬ 
ticularly stimulated by major advances in our ability 
to identify star forming galaxies at high redshifts with 
optical colour techniques (see e.g. Pettini et al 1998 
and references therein). Early indications from com¬ 
parisons of such populations in the Hubble Deep Field 
with similar populations at lower redshift suggested 
that there may have been a peak in star formation 
activity at a redshift of a ~ 2 (Madau et al 1996). 

The theoretical interest in such measurements is 


high, partly because the integrated star formation is 
the major uncertainty in evolutionary models used to 
explain (amongst other things) the local chemical en¬ 
richment in our own galaxy, but also because models 
of galaxy formation based on e.g. hierarchically clus¬ 
tering scenarios can now make specific predictions for 
the integrated star formation rate at different epochs 
(e.g. Pei & Fall, Baugh et al 1997). 

A variety of techniques have been employed to es¬ 
timate the star formation rate at a number of epochs. 
The common feature of such techniques is to iden¬ 
tify a tracer of star formation activity and integrate 
this over volume, correcting for contributions missed 
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through sensitivity limits, obscuration or other incom¬ 
pleteness effects. Tracers that have been employed, 
include: the B-band-luminosity (Lilly et al.) U-band 
luminosity (Madau et al 1996), the Ha luminosity 
(Kennicutt 1983), the mid-IR luminosity (Rowan- 
Robinson et al 1997). The optical tracers are poten¬ 
tially significantly affected by dust-obscuration and 
the corrections for this are uncertain and may be red- 
shift dependent. In contrast Mid and Far-IR luminos¬ 
ity traces the star formation rate in obscured regions 
only. The obscured fraction itself is hotly debated. For 
example, de-reddening of optical-UV spectra does not 
guarantee an unbiased star formation rate, since the 
most heavily obscured and potentially dominant star 
forming regions may be wholly undetected in the UV. 

In this paper we present the first estimate of the 
co-moving star formation density using as a tracer the 
radio luminosity of star forming galaxies. The radio 
luminosity traces the supernovae associated with star 
formation regions, and is not significantly affected by 
dust obscuration. 

The paper is organised as follows. In the first 
section we describe our choice of radio sample. In 
the second we describe our estimation of the radio¬ 
luminosity density including corrections for incom¬ 
pleteness. In the third section we describe the con¬ 
version from luminosity density into star formation 
rate. In the final section we discuss our estimation of 
the star formation rate in the context of other deter¬ 
minations and other implications. 


2 THE 1.4GHZ SAMPLE 

Deep 1.4 GHz radio source counts obtained in the last 
decade (Condon & Mitchell 1984; Windhorst 1984; 
Windhorst et al. 1985) show a steepening below a few 
mjy. This excess of faint sources relative to the Eu¬ 
clidean value can neither be attributed to the giant el¬ 
lipticals and quasars responsible of most of the counts 
observed at fluxes > 5 mJy, nor to the non-evolving 
population of normal spiral and Seyfert galaxies. Op¬ 
tical identification works published so far have shown 
that the sub-mjy sources are mainly identified with 
faint blue galaxies (Kron, Koo & Windhorst 1985; 
Thuan & Condon 1987), often showing peculiar opti¬ 
cal morphologies indicative of interaction and merg¬ 
ing phenomena and spectra similar to those of the 
star-forming galaxies detected by IRAS (Benn et al. 
1993, B93). Although all these works are based on 
very small percentages of identification (due to the 
fact that the majority of the faint radio sources have 
very faint optical counterparts), they all agree that 
most of the sub-mJy radio sources with identifica¬ 
tions brighter than B ~ 22 — 22.5 are star-forming 
galaxies. In contrast, recent extensions to fainter op¬ 
tical magnitudes show strong hints of an increase in 
the fraction of early-type galaxies among the iden¬ 
tifications of sub-mJy sources ( 5 i ,4 ghz > 0.2 mJy; 



Figure 1. The areal coverage 12(5) of the B93 sample as 
a function of 1.4 GHz flux limit 


Gruppioni et al. 1998, in preparation). Thus, down to 
optical magnitudes of B ~ 22.5 most (or eventually 
all) of the star-burst counterparts of sub-mJy sources 
brighter than Si .4 ghz ~ 0.1 — 0.2 mJy should be de¬ 
tected. 

The B93 identification sample consists of 112 
candidate optical counterparts of ~ 500 sub-mJy 
radio sources drawn from the three largest deep 
1.4 GHz radio surveys (0852-1-17 field, Condon & 
Mitchell 1984; 1300-1-30 field, Mitchell & Condon 
1985; 0846-1-45 field, Oort 1987). All the selected 
sources have radio fluxes greater than 0.1 mJy and op¬ 
tical counterparts brighter than V = 20.0 (0846+45 
field), B = 21.4 (0852+17 field) and B = 22.3 
(1300+30 field). The radio flux limit varies across the 
survey and the areal coverage {Q{S)) of the B93 sam¬ 
ple is a complicated function of radio flux, which is 
plotted in Figure Spectra have been obtained for 
87 of these, providing object type and spectroscopic 
redshifts, though at 2 ; > 0.35 Ha is unobservable in 
their spectra so the spectroscopy is only complete at 
low redshifts. Of these sources with spectra, 47 turned 
out to be star-burst galaxies. The B93 sample is the 
largest spectroscopic sample of sub-mJy sources so 
far available in literature and, due to its radio and 
optical limits, is the most complete sample of radio- 
selected sub-mJy star-forming galaxies. This sample 
is therefore the best currently available and is the 
most suitable for estimating the star-formation his¬ 
tory as traced by the radio luminosity density of star 
forming galaxies. 

In all the following analysis we use only those ob¬ 
jects within B93 which Rowan-Robinson et al. iden¬ 
tified as star-burst galaxies. 


3 THE LUMINOSITY FUNCTION AND 
LUMINOSITY DENSITY 
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3.1 Basic Method 

To calculate both the luminosity function and the lu¬ 
minosity density we require an estimate of the effec¬ 
tive volume within which each object of the sample 
could have been observed. Besides a « < 0.35 spec¬ 
troscopic redshift limit, the B93 sample has two flux 
limits, one radio and one optical, this effective volume 
is thus given by the following expression 

Vi = J P{0\Li,LB,uz)n^^^—dz ( 1 ) 

where P{0\L, Lb, z) is the probability of observing 
an object of radio luminosity L and optical luminosity 
Lb Sit redshift « within our survey, given by 

= VM^^P(C9|Ss), (2) 

^ ^max 

where the summation is over each of the three sur¬ 
vey areas, in our case P{0\Sb) is either 1 or 0 
depending on whether the magnitude is below or 
above the limit in that survey area. The radio fluxes 
are calculated assuming a radio spectral index, a = 
cilog(S',/)/dlog(i/) = 0.8; the optical fluxes are calcu¬ 
lated using the if-corrections for Scd galaxies given 
by Metcalfe et al. (1991) and we assume a cosmol¬ 
ogy of flo = 1, Ha = 0, Po = SO/iso km s“^ Mpc“^ 
throughout. 

Having determined Vi we can use the estimator 
tp{L) = 'Y^'i/Vi (Schmidt 1968), and likewise esti¬ 
mate the luminosity density by 'l'i ,4 = '^LijVi. In 
the absence of luminosity cuts these estimators are 
unbiased, however one might worry about using these 
estimators to determine the radio luminosity func¬ 
tion and luminosity density if the optical flux cut was 
much more severe than the radio, leading to poor sam¬ 
pling of the radio luminosity plane. We will demon¬ 
strate that this is not the case. 

In the presence of luminosity cuts, introduced for 
example by redshift cuts these estimators need to be 
modified thus: 


■I?!.4 


C{L 


min^ 



(3) 


where C{Lmi-n, Lb, min) is a measure of the “com¬ 
pleteness” of the sample being the fraction of lumi¬ 
nosity density above the minimum observable radio 
and optical luminosities i.e. 


f+ao r+oo Lipi^L^Lg-^dlogLBdlogL 

/_^r /_"T^^(^--^s)dlogLsdlogL 


tp{L,LB) is the 1.4 GHz optical bi-variate luminosity 
function. 

We will show that these correction factors are 
principally determined by the shape of the lumi¬ 
nosity function. Maximum likelihood estimators (e.g 
Sandage et al 1979) can determine the shape of the lu¬ 
minosity function independently of density variations 
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Figure 2. Optical Luminosity versus Radio Luminosity 
for the Condon sample (crosses) and for the B93 sample 
(filled circles) and for the B93 sample corrected to z = 0 
luminosity according to L = L{z = 0)(1 + (filled 

squares). Over-plotted is the locus of the radio and optical 
flux limits for different redshifts, no evolution upper, de- 
evolved to z = 0, lower 


which can strongly effect the XjV estimator. Using a 
method akin to the maximum likelihood technique we 
take some reasonable a priori models for '^{L) (and 
selection effects given by Pi{0\L)) to predict the ex¬ 
pected luminosity distribution independently of den¬ 
sity fluctuations, which can be compared with the ob¬ 
served luminosity distribution using a K-S test. The 
expected luminosity distribution for an object in the 
survey is given by 


f(L) = P^{0\L)9{L) 

Pi{0\L)^(L)dlogL 


(5) 


The expected luminosity distribution of the sample is 
then given by 


3.2 Impact of optical flux limit 

To demonstrate that the optical selection effects are 
less significant than the radio selection effects in Fig¬ 
ure 1^ we plot the radio and optical luminosities of the 
Condon (1987, C87) spirals, which has > 99% com¬ 
plete radio detections in a sample of Revised Shapley- 
Ames galaxies (see Condon 1987 for more details). 
Note that there is effectively only one flux limit, so 
this correlation is not a “distance vs. distance” selec¬ 
tion effect. 

The z < 0.35 star-bursts from B93 are also plot- 
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ted in Figure as is the passage of its double flux 
limits to 2 = 0.35. 

The immediate point to notice from this Figure 
is that the L, Lb plane is not well sampled by the 
B93 sample, the low radio luminosity region is poorly 
sampled. However this is due to the small volume at 
low redshift, there appears to be reasonable sampling 
across Lb for a give radio luminosity. 

For a non-evolving population, virtually all the 
radio sources lie above our optical limit at 2 ; < 0.35. 
In other words, if we regard our sample as a purely 
radio flux limited sample, then we can be confident of 
very high completeness to the radio flux limit. Rowan- 
Robinson et al. demonstrated that luminosity evolu¬ 
tion of the same strength as required for IRAS galax¬ 
ies L{z) = L{z = 0)(1 -I- z)^'^ was sufficient to ex¬ 
plain the sub-mjy radio counts. Such a pure luminos¬ 
ity evolution is equivalent to a large negative radio K- 
correction term, and the corrected positions in Figure 
^ also demonstrate a high completeness. 

We can quantify this by using the radio-optical 
correlation to estimate the incompleteness as a func¬ 
tion of redshift, then use the observed redshift dis¬ 
tribution to get a first-order estimate of the overall 
incompleteness: 

/ = 100% X iVobs/ A:(Ln,in,i, Zi) = 98.4% (6) 

where / is the completeness. Nobs the size of the 
star-burst sample. A: is 1 over the fraction of the 
radio-optical correlation above the optical flux limit 
as a function of redshift and luminosity, imin.i is the 
faintest radio luminosity observable at redshift Zi, and 
the sum is performed over the star-burst sample. The 
B93 sample does appear to have a slightly larger dis¬ 
persion in the radio-optical plane than the local sam¬ 
ple, but this larger scatter still yields 96.7% complete¬ 
ness. We conservatively assume no passive stellar evo¬ 
lution, the effect of which would be to increase our 
completeness further still. N.B. this measure of “com¬ 
pleteness” is not the same as that described in Equa¬ 
tion The assertion that the optical selection does 
not have a significant additional impact on the radio 
selection is consistant with the findings from deep op¬ 
tical observations of the Marano Field (Gruppioni et 
al. in preparation) that the faintest optical identifica¬ 
tions of sub-mJy radio sources do not appear to be 
star-burst galaxies. 


3.3 1/y Luminosity Function 

The luminosity function is plotted in Figure which 
agrees with the calculation in Rowan-Robinson et al. 
(1993). It is clear from Figure ^ that the faint end 
of the luminosity function is not well-constrained by 
this sample. This is due to the small volume at low 
redshift of this sample. For 2 < 0.05 (where the limit¬ 
ing radio luminosity is L = 1.1 x 10^^h^Q^WHz“^ the 
B93 sample covers only 2700h^Q Mpc^, and using the 


Starburst LF 



log.oPi.4 (W/Hz) 


Figure 3. The 1/F Luminosity function for the B93 sam¬ 
ple: filled circles. Also illustrated is the luminosity function 
from Condon 1987, upper curve, the warm IRAS luminos¬ 
ity function of S90, lower curve. Also illustrated is the 
estimate of luminosity function from Windhorst (1985) 


S90 “warm” luminosity function we would only expect 
1.75 galaxies within this volume. This small volume 
at low redshifts is also very susceptible to large scale 
structure variations. Assuming the power spectrum 
of Peacock and Dodds (1994) we would expect RMS 
density fluctuations of around 60 to 60 per cent over 
such a volume (allowing for the survey being spilt into 
three). 

The co-moving number densities are consistent 
with the single star-burst at 2 < 0.35 detected in 
the Hubble Deep Field (Richards et al. 1998), as well 
as with the luminosity function from the 2 < 0.35 
Marano star-bursts (Gruppioni et al 1998, in prepa¬ 
ration). This data is however in strong disagreement 
with the very high star-burst density claimed in the 
micro-Jansky radio population by Windhorst et al. 
(1995). Figure ^ plots the 2 < 0.35 Windhorst et al. 
star-bursts, which have an order of magnitude higher 
space density than our estimate. The nature of the 
fijy population is controversial however, with Ham¬ 
mer et al. (1995) claiming the population is domi¬ 
nated by weak AGN, albeit in a sky area which may 
be biased by the presence of a cluster. Part of the dis¬ 
crepancy may also be due to high frequency (5 — 8.5 
GHz) selection in both the Hammer and Windhorst 
samples, where the radio spectral index may not be 
steep due to the flat-spectrum contributions from Hll 
regions and possibly AGN. The high-frequency se¬ 
lected sub-mjy//ijy population may therefore be be¬ 
ing sampled higher up the luminosity function, or may 

© 1998 RAS, MNRAS 999, 0-0 





5 


be contaminated by Seyfert galaxies. This justifies our 
choice of a low-frequency selected sub-mjy radio sam¬ 
ple, albeit in retrospect. 


3.4 Testing the shape of the luminosity 
function 

We have already discussed that the B93 sample is in¬ 
sufficient to determine ip{L) at faint radio luminosities 
since the volume and thus number of objects at low 
redshift is small. 

To estimate the corrections to our estimate of the 
luminosity density (Equation^ we are only interested 
in the shape of the luminosity function (where shape 
includes a measure of the typical Luminosity e.g. L*) 
and not the overall normalisation which cancels in 
Equations ^ and p. We thus wish to test the shape 
of B93 luminosity function against some reasonable a 
priori functions. 

The IRAS luminosity function is well determined 
even at the faint end, it is natural to test an IRAS 
luminosity function modihed via the FIR/radio cor¬ 
relation (Helou, Soifer & Rowan-Robinson 1985). For 
a galaxy with the model star-burst spectrum of Efs- 
tathiou, Rowan-Robinson & Seibenmorgen (1998, in 
preparation) this relation can be translated to 


Sqo ^rn 

= 133 5i,4 GHz 

(7) 

or 

Lqo jjm 

= 133 Li,4 GHz 

(8) 


We consider two models, the first (11) is the 
“warm” luminosity function of S90(their solution 25), 
(modified so L* = 1331/60*). It is appropriate to use 
the “warm” sample, since these are the actively star¬ 
forming galaxies, similar to this sample. This luminos¬ 
ity function is also clearly similar in shape and nor¬ 
malisation to our observed radio luminosity function 
(Figure 1^. The use of this model luminosity function 
is independent of the overall normalisation, and as¬ 
sumes the shape is constant with redshift. It is known 
that IRAS galaxies strongly evolve (e.g. S90, Oliver et 
al 1994, etc.), and if this evolution only changes the 
normalisation of the luminosity function then this is 
immaterial. However if shape of the luminosity func¬ 
tion was to change with redshift, this would affect 
our estimate of the luminosity density. Neither IRAS 
samples, nor this sample alone are sufficient to deter¬ 
mine what sort of evolution is appropriate to these 
types of samples; however the faint radio counts do 
not support pure density evolution as strong as that 
suggested in S90 (Rowan-Robinson et al. 1993). Thus 
it is appropriate to consider some model with an evo¬ 
lution of the Luminosity with redshift. Our second 
model 12 assumes the same zero redshift luminosity 
function as Z1 but has L = Lo(l + ^)^'^. Since the S90 
2 = 0 luminosity function was not calculated with 
this evolution, this model is not entirely consistent 
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and will over-predict the numbers of luminous galax¬ 
ies, and lead to an under-correction to the Luminosity 
density, nevertheless the two models together give a 
reasonable idea of the level of uncertainty in the lu¬ 
minosity function, and would bracket a more natural 
model where Luminosity evolution was included in 
the estimate of the Luminosity function self consis¬ 
tently. 

Our third model, 13, assumes the shape of the 
low 2 Condon et al (1987) radio luminosity function, 
with no evolution. 

In addition to the luminosity functions we con¬ 
sider a number of different models for the effect of the 
joint optical and radio limits on the observable radio 
luminosities Pi{0\L)-. 


Pi{0\L) 
PiiO\L) = 


p{0\L) = 

^max 

^{S{L, Zi)) flB{SB{LLB,i/Li, Zi)) 
^max ^S,max 

^{S{L,Zi)) nBPiLB\L)dlgLB 


n 


max 




,max 


(9) 

( 10 ) 

( 11 ) 


(In the above equations B represents the optical 
which may be either B or V band, where correspond¬ 
ing magnitudes of individual objects are computed 
using the S5000 /•S'tooo ratio in B93 and magnitude lim¬ 
its are adjusted using the average ms — mv = 0.7.) 
Model mO (Equation 0) assumes that the optical mag¬ 
nitude of the object is totally unrelated to the radio 
luminosity so the only selection that comes into play 
is the radio flux limit. If we knew the radio flux limit 
as a function of position Equation ^ would reduce to 
Pi{0\L) = 0, L < Z/inin,i; Ei(Cl|Z/) = 1, L > L/min.i. 
Equation ^ (model ml), assumes that the lower lu¬ 
minosity objects potentially missed through the opti¬ 
cal selection have the same radio/optical colour as the 
object under consideration, this will be satisfactory if 
the colours within the sample are representative of 
the colours in the underlying population (in fact the 
optical and radio flux limits will survey to restrict 
the range of colours); in our case of simple magnitude 

limits — 0, ,Sb Bb. liin, ^^B/^B,inax — 

1, Bb > S'B.iim. Model m2 (and m3) take the form 
of Equation ^ and require some estimate for the 
marginal optical radio distribution function P(Lb|L): 
model m2 takes this from Saunders et al. (1990, S90), 
assuming Leo = 1331/; model m3 takes this from the 
Condon et al (1991). 

Using the different models (mO, ml, m2, m3) for 
the optical and radio selection criteria Pi{0\L) and 
the three models for the luminosity function 11,12,13 
we calculate the luminosity distributions and compare 
with the observed luminosity distributions using a K- 
S test, these results are summarised in Table [ij (In 
this and subsequent analysis we have excluded the 5 
objects which fall below the optical completeness lim¬ 
its quoted in Benn et al..) From this we can see that 
none of the models tested can be significantly ruled 
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Pi(0\L) 

mO 

ml 

m2 

m3 

ip{L) 

n 

11 

11 

11 

P(KS) 

0.13 

0.64 

0.13 

0.13 

Pi{0\L) 

mO 

ml 

m2 

m3 

ifi(L) 

12 

12 

12 

12 

P{KS) 

0.30 

0.83 

0.29 

0.30 

Pi(0\L) 

mO 

ml 

m2 

m3 

ip{L) 

13 

13 

13 

13 

P(KS) 

0.04 

0.33 

0.04 

0.04 


Table 1. Comparison of the observed and expected lumi¬ 
nosity distributions under various assumptions about the 
underlying luminosity function and selection effects. 



Logic (Li,i GHz/1^'*' Hz ‘0 


Figure 4. The observed and expected luminosity distribu¬ 
tions. Solid histogram, observed distribution. Dotted curve 
expected luminosity distribution assuming “warm” IRAS 
luminosity function (11) and no accounting for optical se¬ 
lection (mO, similar to m2, m3). Solid curve ll and op¬ 
tical selection determined by objects colours (ml). Dot- 
dashed curve, evolving “warm” IRAS luminosity function 
(12). Dashed curve Condon Luminosity function (13). 


out. There is very little difference between the mod¬ 
els which assume either no optical selection effects or 
optical selection effects determined from a priori esti¬ 
mates of the radio/optical distribution, this is because 
these estimates do not predict any significant loss of 
sources to the optical selection effects. A better fit is 
provide by models which assume the optical selection 
is determined by the distribution of colours within the 
sample itself (ml). All forms of luminosity function 
give acceptable fits, assuming ml. 


3.5 The luminosity density 

The correction factor {C{L min, Lb, min), Equation ^ 
to be applied to the luminosity density in the presence 
of luminosity (i.e. redshift) cuts, can be re-expressed 
and estimated in a number of ways, we choose to ex¬ 


press it as 

/+“ Lip{L)dlogL 

^ ^min _ 

/+^L(p(L)dlogL 

r + oo r + oo l^g d log L 

/10^ 

L(p(L, LB)d log Lsdlog L 

The last term in this equation represents the loss 
in luminosity density due to the optical cuts alone. We 
have calculated this term at 2 = 0.05 and 2 = 0.175 
using the Condon bivariate luminosity function, the 
Saunders bivariate luminosity function, and for a ver¬ 
sion of the Condon bivariate luminosity function with 
increased dispersion appropriate for the B93 sample. 
The largest value we get for the second term is 010“®, 
which is negligible compared to the former term in¬ 
volving only the radio luminosity density Ol. Since 
the normalisation cancels out, the only significant cor¬ 
rection to be applied is due to the shape of the lumi¬ 
nosity function. 

With surveys of small volume density fluctua¬ 
tions on the scale of the survey volume can contribute 
significantly to the uncertainty in derived quantities. 
To estimate the errors due to this cosmic variance 
we determine the effective wavenumber k appropriate 
for a cube with volume one third of the total volume 
in the redshift range, since the survey is composed 
of three independent volumes. From this we can es¬ 
timate the density fluctuations — A^(fc) from the 
power spectrum of Peacock and Dodds (1994), these 
we reduce by a factor of three, assuming the three 
separate areas are independent. These errors are ap¬ 
proximate and do not take properly into account the 
survey geometry or the fact that lower luminosity ob¬ 
jects sample smaller volumes, all these effects would 
serve to increase the errors. These errors are slightly 
larger than the shot noise terms to which they are 
added in quadrature. We tabulate the cosmic variance 
errors for this and other samples used for estimating 
luminosity densities(Table |^). 

In Table ^ we tabulate the un-corrected luminos¬ 
ity density for the B93 sample over 4 redshift ranges. 
We also tabulate the correction factors that need to be 
applied assuming various luminosity functions which 
we have shown to be compatible with this sample, (for 
the evolving luminosity function we have take then lu¬ 
minosity function shape to be defined at roughly the 
centre of the redshift bin). We also tabulate the lu¬ 
minosity density estimated from a Schecter function 
fit to the luminosity function Condon RSA excluding 
AGN (Serjeant et al, in preparation). We also tabu¬ 
late the luminosity density for warm IRAS galaxies 
from S90, scaled by "heoMm = 133'l>i.4. 

We can see from this that the total luminosity 
density in the B93 sample as a whole is in good agree¬ 
ment with the luminosity density from the Condon 
sample. The luminosity density in the higher redshift 
bin is larger than that in the lower redshift bin by a 
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Sample 

Area 

sub. 

Zmini Zmax 

Ka/h-so 

a 

z 

4 (Ha) 


/ 10 -'‘sr 



/Mpc-i 



/h5oWMpc-3 

B93: 

3.4 

3 

0, 0.05 

0.30 

0.72 

0.0-0.35 

31.8627 

B93: 

3.4 

3 

0.05, 0.175 

0.095 

0.41 

0.05-0.35 

31.7764 

B93: 

3.4 

3 

0.175, 0.35 

0.055 

0.26 

0.05-0.175 

31.7554 

B93: 

3.4 

3 

0.05, 0.35 

0.052 

0.25 

0.175-0.35 

31.2052 

B93: 

3.4 

3 

0, 0.35 

0.052 

0.25 



Tr 0 ss 6 & 

0.42 

5 

0, 0.3 

0.040 

0.16 

Table 4. Ha Luminosity density estimated from star- 

A/I&dciox! 






forming galaxies within the B93 sample. Errors and cor¬ 







rection factors will be similar to those in .3 

Lilly 

0.42 

5 

0.2, 0.5 

0.028 

0.10 



Lilly 

0.42 

5 

0.5, 0.75 

0.024 

0.088 



Lilly 

0.42 

5 

0.75, 1 

0.023 

0.079 

included the spectroscopic limits in our effective vol- 

Gallego 

1400 

1 

0, 0.045 

0.0089 

0.036 

time. 


Gronwall 

210 

1 

0, 0.085 

0.0092 

0.017 




Table 2. Cosmic variance errors (tr = ^ A^(fc)) within 
surveys of different areas, redshift ranges and different 
number of subdivisions. For the B93, Tresse & Maddox 
and the low z Lilly samples these errors are comparable 
with the shot noise errors. 


z 

log 41,4 

log 

Rmin) 


/ft,5oWHz-iMpc-3 

11 

12 

/3 

0-0.05 

19.27 ±0.14 

0.0 

0.0 

0.0 

0-0.35 

19.09 ±0.14 

0.0 

0.0 

0.0 

0.05-0.35 

19.09 ±0.14 

0.03 

0.02 

0.05 

0.05-0.175 

18.91 ±0.20 

0.03 

0.02 

0.05 

0.175-0.35 

19.01 ±0.16 

0.21 

0.12 

0.42 

0 - 0.1 

18.69 ±0.05 

0.0 

0.0 

0.0 

0 - 0.1 

19.08 ±0.04 

0.0 

0.0 

0.0 


Table 3. Estimates of the raw 1.4GHz luminosity density 
in different redshifts ranges and corrections that need to be 
applied under different assumptions about the underlying 
luminosity function. The 0 < z < 0.05 luminosity density 
is taken from Serjeant et al. (in preparation). The last two 
rows are obtained from converting the “warm” and total 
IRAS luminosity densities of S90 to a 1.4Ghz luminosity 
density using 133Z/i,4 = Leo- 


factor of between 1.7 and 3.7, depending on the lu¬ 
minosity function applied, consistent with luminosity 
evolution at a rate of (1 -|- however, due to the 

cosmic variance this is not very significant. 

The luminosity density from the low-redshift 
sample alone is a factor of around two lower than the 
Condon sample and similar to the luminosity density 
as estimated from the “warm” IRAS galaxies. Due 
to the cosmic variance the discrepancy with the Con¬ 
don sample is not signihcant, while the discrepancy 
between the Condon sample and the “warm” IRAS 
sample is. It is likely that the Condon sample also in¬ 
cludes the radio equivalent of the IRAS “cool” Cirrus 
galaxies. 

Using the Ha luminosities quoted in Benn et al 
1993, corrected for slit loss using the ARM magni¬ 
tudes and continuum fluxes we can use the same ef¬ 
fective volumes to estimate the Ha luminosity density 
within this sample. Table ^ This may be an underes¬ 
timate of the total Ha luminosity density since from 
some objects Ha is not detected and we have not 


4 CONVERSION FROM LUMINOSITY 
DENSITY TO STAR FORMATION 
RATE 


Nearly all of the radio emission from star-forming 
galaxies is synchrotron radiation from relativistic elec¬ 
trons and free-free emission from Hll regions. Only 
massive stars (i.e. M > 5Mq) ionise the Hll regions 
and produce supernovae, whose remnants accelerate 
most of the relativistic electrons. Such massive stars 
have lifetimes much shorter than the Hubble time, 
so the current radio luminosity is proportional to the 
recent star formation rate (Condon 1992). The su¬ 
pernova rate is directly related to the non-thermal 
radio luminosity, which at 1.4 GHz is dominant over 
the thermal component (see Condon & Yin 1990 for 
the most reliable derivation of this relation, calibrated 
with Galactic supernova remnants). Moreover, since 
all stars more massive than M = 5Mq become radio 
supernovae, the radio supernova rate is determined 
directly by the star formation rate. Thus, the star 
formation rate can be obtained by non-thermal radio 
luminosity, following Condon (1992): 


SFR(M > 5 M 0 ) 


L1.4 GHz [WHz-^] 
5.3xl02i(^)- 


M^yr ^(13) 


where u is the frequency and a is the non-thermal 

radio spectral index (as defined above). If we assume a 

Salpeter initial mass function (IMF; tp{M) oc 

0.1 — 125 Mq) we obtain for the total star formation 

rate 


SFR(M > O.IM 0 ) = 


Li.4 gh. [WHz-^] 

7.63 X 1020 


(14) 


In Figures & ^ we convert our radio luminosity 
densities (Tabled, row 3, correction 11) to star forma¬ 
tion rates together with that from the RSA sample 
(Table P, row 1). 

It is instructive to compare with the U-band lu¬ 
minosity densities at redshifts z i 1, from Treyer et 
al. (1998) and the Canada-France Redshift Survey 
(CFRS, Lilly et al. 1996). For the ultraviolet conver¬ 
sions we assume 


SFR{M > O.IM 0 ) 


L = [Wl 
2000A '■ ‘ 

7.94 X 1020 


(15) 
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0.0 0.2 0.4 0.6 0.8 1.0 

Redshift 



0.0 0.2 0.4 0.6 0.8 1.0 

Redshift 


Figure 5. The co-moving star-formation rate 2 < 1 as es¬ 
timated from: Filled circles = 1.4GHz (Left to right: Con¬ 
don, B93) Open circles = Ha (Gallego, Gromwell, B93, 
Tresse &: Maddox) Open triangles = UV (Treyer, Lilly, 
Lilly, Lilly). Corrections for reddening have not been made 
or have been removed 

L • [W] 

SFR{M > O.lMo) = (16) 

taken from Treyer et al. (1998) and Cowie et al. (1997) 
respectively. In Figure la we do not apply a reddening 
correction, in Figure ^ the luminosity densities are 
de-reddened following Treyer et al. (1998). 

Also plotted are the Ha luminosity densities from 
the CFRS (Tresse & Maddox 1997) the KISS sample 
(Gromwell 1998) and from the local {z = 0) sample of 
Gallego et al. (1995), converted to total star formation 
rates for this IMF using 

5Ti?(M>O.lM0) = ^^l^ (17) 

from Madau et al. (1996). Both Ha luminosity densi¬ 
ties were corrected for reddening using Balmer decre¬ 
ments and assuming a simple dust screen, in Figure ^ 
we remove the average reddening correction ~ 1 mag 
(Tresse & Maddox 1997). 

A priori we would have expected the 1.4GHz to 
provide a higher estimate of the star-formation rate 
than an un-reddened optical estimate. This is be¬ 
cause the radio should be measuring the total star- 
formation rate which the optical only measures the 
un-obscured rate. At the low redshift (0 < 2 < 0.05) 
this does indeed appear to be the case, the Gondon 
luminosity density is well above that estimated from 
the Ha of Gallego et al. (1995) and Gromwell (1998), 
and above the UV estimate of Treyer et al. (1998). 
At intermediate redshifts (0.05 < 2 < 0.35), probed 
by the B93 survey, the Tress & Maddox Ha estimate 
from the CFRS and the UV estimate of Treyer et al. 
(1998), the estimates agree, if the reddening correc- 


Figure 6 . The co-moving star-formation rate 2 < 1 as es¬ 
timated from: Filled circles = 1.4GHz (Left to right: Con¬ 
don, B93) Open circles = Ha (Gallego, Gromwell, B93, 
Tresse & Maddox) Open triangles = UV (Treyer, Lilly, 
Lilly, Lilly) Corrections for reddening have been made or 
retained 

tions are applied then the optical estimators can even 
exceed the radio estimate (Figure^. This is initially 
surprising until we notice that the errors are large, we 
note that we have applied the cosmic variance terms 
to all the data points and these are significant, partic¬ 
ularly for the B93 sample and the Tresse & Maddox 
sample. We thus suspect that the coincidence of these 
various estimators may be due to differences in the 
mean density of these survey volumes which masks a 
difference in their sampling of the star-formation rate. 
This assertion is strengthened when we see that the 
Ha estimate of the star-formation rate from the B93 
sample itself is further below the radio estimate than 
the Tresse & Maddox point is above it. 

An additional star-formation indicator is the far 
infrared. Several authors (Scoville & Young 1983; 
Thronson & Telesco 1986; Condon 1992; Rowan- 
Robinson et al. 1997) have derived relations to convert 
from 60 fim luminosity to star formation rate. These 
calculations assume some proportion of the optical 
luminosity produced by young stars is absorbed and 
re-emitted in the infrared. (This same extinction fac¬ 
tor should naturally be applied to the star-formation 
rates calculated from the UV.) From such calcula¬ 
tions it is of course possible to deduce the FIR/Radio 
correlation. Condon (1992) deduced the FIR/Radio 
correlation assuming an obscuration factor of 2/3., 
recently Cram et al. (1998) seems to indicate a close 
agreement between the radio SFR estimate and the 
far-infrared one obtained considering ~ 100 per cent 
reprocessing of starlight in the far-infrared. In general 
it does appear that the extinction fraction is high. It 
would in fact be possible to use the FIR/Radio cor- 
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relation to determine the extinction fraction for any 
assnmed underlying IMF, this could then be feed self 
consistently into the UV and Ha SFR estimates. Such 
an analysis is beyond the scope of this paper. 

Instead we use the FIR/Radio correlation sim¬ 
ply to convert local estimates of the 60pim luminosity 
density from S90 to 1.4GHz luminosity density and 
then use the 1.4GHz calibration above which is in¬ 
dependent of extinction. We plot in Figure the 
star-formation rate estimated thus from the S90 lu¬ 
minosity densities for both “warm” IRAS galaxies and 
for all IRAS galaxies. It is not clear that the “cool” 
IRAS galaxies trace star-formation as do the warm 
IRAS galaxies, since the cool emission is from cir¬ 
rus clouds which may be illuminated by older stellar 
populations (similarly the UV density of quiescent or 
weakly star-forming galaxies might also be dominated 
by old stars). 


5 CONCLUSION 

We have estimated the 1.4 GHz luminosity density 
from a sample of radio star-burst galaxies z < 0.35, 
demonstrating that the optical selection criteria ap¬ 
plied to this sample are not important to this deter¬ 
mination. The small volume of this survey ai z < 0.05 
prevents us from having accurate estimates of the 
luminosity density from low radio luminosity star¬ 
forming galaxies. Nevertheless we demonstrate that 
a number of reasonable a priori luminosity functions 
can be used to estimate the luminosity density missed. 
These corrections are small for the sample as a whole, 
though larger and more disparate for a higher redshift 
sub-sample. The 1.4 GHz luminosity density increases 
as we move to higher redshifts, though the strength 
of this effect is dependent on the assumed luminos¬ 
ity function and may not be conclusive with a sample 
of this volume. Overall the 1.4 GHz luminosity den¬ 
sity agrees with that obtained from 1.4 GHz measure¬ 
ments of the RSA sample of Condon et al.. 

The 1.4GHz estimator of star-formation rates 
should be unaffected by extinction, whereas local es¬ 
timates from the UV and Ha should be. We were 
thus surprised to find the estimate from the radio 
broadly consistent with these other estimates, par¬ 
ticularly since the local estimate from the Gondon 
luminosity density was considerably higher than from 
optical measures. Since the luminosity density esti¬ 
mated from Ha within the B93 sample produces a 
much lower estimate of the star-formation rate than 
the 1.4GHz estimate and the (optimistic) cosmic vari¬ 
ance errors in the B93 and CFRS samples are com¬ 
parable to the shot noise errors we attribute this ap¬ 
parent agreement to clustering fluctuations. We thus 
strongly suggest that spectroscopic follow-up of larger 
area radio surveys is needed to accurately determine 
the un-obscured star-formation rates 0 < z < 0.3. 

The star formation rates at these redshifts do not 


conflict with the Far-IR background measurements, 
since the sub-mjy radio population make a negli¬ 
gible contribution to the FIR background. At these 
redshifts the 1.4 GHz and 850/im fluxes are roughly 
equal, and the Puget et al. (1996) 850/rm background 
corresponds to about 26 Jy per square degree. 
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